Anomalous fluctuations of active polar filaments 
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Using a simple model, we study the fluctuating dynamics of inextensible, semiflexible polar 
filaments interacting with active and directed force generating centres such as molecular motors. 
Taking into account the fact that the activity occurs on time-scales comparable to the filament 
relaxation time, we obtain some unexpected differences between both the steady-state and dynamical 
behaviour of active as compared to passive filaments. For the statics, the filaments have a novel 
length-scale dependent rigidity. Dynamically, we find strongly enhanced anomalous diffusion. 



Filamentous proteins are major components of the cell 
cytoskeleton JD. Examples are actin filaments, micro- 
tubules and intermediate filaments. Their mechanical 
properties are important for cell stability and support 
and have been well studied at thermodynamic equilib- 
rium by in vitro experiments. However, the conditions 
in the living cell are very different from those in the lab- 
oratory. Protein filaments interact with other proteins 
such as molecular motors and cross-linkers. This has led 
to a flurry of recent in vitro experiments of mixtures of 
filaments and their associated proteins in order to com- 
pare with their purified state |§-0). From a theoretical 
point of view, the proteins are typically far from equilib- 
rium and therefore even to understand their steady-state 
behaviour, one has to study their dynamics. Non equilib- 
rium effects have also recently been studied in biological 
membranes ||. 

In this paper, motivated by recent experiments on F- 
actin and myosin in the presence of ATP [El, we study 
one example of the non-equilibrium behaviour of bio- 
filaments: the fluctuating dynamics of polar filaments 
with active centres. This is also a model system for the 
study of non-trivial aspects of semiflexible filament dy- 
namics. A key point of our analysis is the fact the activity 
of the proteins occurs over a time-scale r which may be 
comparable to the relaxation time of the filament. Un- 
like recent work on motile solutions B , the active centres 
considered here are associated with single filaments (i.e. 
are not cross-links) and cannot move one filament with 
respect to another. The viscosity of the solvent is given 
by 77. We note that the typical energy scale of a biochem- 
ical reaction is of the order of a few ksT at physiological 
conditions. 

The filament can be parameterised by a curve through 
its centre, R(s) (sec Fig. |l|). The unit tangent- vector 
is defined as t(s) = dH/ds. An easily measured quan- 
tity using e.g., video microscopy is the steady-state tan- 



gent correlation function Ct(s) = (t(s) ■ t(0)). For a 
semiflexible polymer in equilibrium it decays exponen- 
tially, Ct(s) = exp [— |s|/i p ] which defines the persis- 
tence length L p . 

Main results We calculate Ct(s) and find that due to 
activity on a time-scale r, the filaments develop a length- 
scale dependent bending rigidity. On short length-scales, 
typical conformations have the bare persistence length 
whilst on longer length-scales the filaments may be char- 
acterised by a lower 'renormalized' persistence length. 
There is a cross- over length between the two regimes, 
£ c ~ {TkBTLp/2-Kr)) 1 ^ . Therefore an analysis of the 
filament conformations can give information about the 
time-scale t at which the activity occurs. We also obtain 
modified relaxational dynamics with anomalous diffusion 
which we show has an effect on the high frequency vis- 
coelasticity. The shear modulus of the active filament 
solution apparently has a frequency dependent effective 
temperature. There is a crossover from at high frequen- 
cies a shear modulus corresponding to the bare tempera- 
ture to at lower frequencies a modulus equivalent to the 
higher renormalized temperature both with a scaling of 
w 3 / 4 . In the cross-over regime, the modulus appears to 
have a stronger frequency dependence with a power law 
G*(u>) cx oj a where a > 3/4. A simple consequence of 
filament polarity is a small ballistic component to the 
motion of the filament, 

Explicit analysis We model the filaments with 
the Kratky-Porod worm-like chain which takes 

into account the bending energy cost of the chain. 
The worm-like chain Hamiltonian is 7i tiI ; c [{R(s)}] = 



§ /-l/2 ^ s (l^) • ^he stretching energy of chain 
molecules is much higher than the bending energy and 
we may consider the chain as inextensible. The persis- 
tence length is defined as L p = n/ksT. 

The dynamics at finite T may be expressed by the 
Langevin equation 



d_ 

dt 



R(s,t) - J ds' H(s,s') 
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V(s,f) = f(M) + f (m) (M), 



(1) 
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where A(s, t) (a Lagrange multiplier) is an instan- 
taneous, fluctuating 'tension' which enforces the local 
inextensibility of the chain. The inextensibility con- 
straint fixes the tension, A(s,t), by = |^| — 1. 
Because of this constraint, the dynamics of semiflexi- 
ble filaments is non-linear and generally insoluble. In 
addition to the thermal velocities f(s,t), there is an 
additional active velocities f( m )(s,i) correlated over a 
time r. The polarity of the filament implies the active 
forces have a mean direction along the contour of the 
filaments (see Fig [l]). Hydrodynamics is taken into ac- 
count by the mobility tensor H(s, s'). In the Stokes ap- 
proximation valid for low Reynolds number, H(s, s') — 
H(R(s) — R(s')) and generates a non-dissipative term ||, 
Vi(s,t) = kgT ds'5Hij(s, s')/5Rj(s' ,t). In addition 
f(s, t) will also depend on R(s,i) ||. We assume that 
the activity does not affect the thermal forces, i.e. the 



active proteins do not affect the collisions of the solvent 
molecules with the filament. Then the thermal velocities 
satisfy the fluctuation-dissipation theorem, and f(s,i) 
have zero mean and Gaussian fluctuations, given Q by 
(f i (s,t)f j (s',f))=2k B TH ij (s,s')8(t-t'). 

We consider rod-like segments and restrict the anal- 
ysis to length-scales I below L p so that we can decom- 
pose the dynamics into transverse and longitudinal mo- 
tion (see Fig. |l|) and write the position of the fila- 
ment as R(s, t) = (s — rii (s, t))u(t) + r±(s, t) where s S 
{— 1/2, 1/2} and u is a time-dependent unit vector giving 
the orientation. We can obtain results within a system- 
atic small gradient expansion for |rj_(s, t)\, |r*|| (s, t)\ <C 
s. The mobility tensor is given by Hij(s,s') — h(s — 
s') {{8ij — UiUj) + 2iiiUj + 0(\d s r±\ 2 )} where in this and 
the following, for a function A(x), d x A = dA/dx. We can 
thus decompose eqn. (Q) into parallel and perpendicular 
components in an expansion to 0(|<9 s rj_| 2 ), 

(2) 
(3) 



d t r x (s,t) = J ds'h(s-s') [- K 5 4 ,r ± + A(s' ,i)d 2 ,r ± + d s ,Ad s >r ± ] + f ± (s, t) + f[ m) (s, t) 
d t r\\{s,t)= f ds^his-s') [-«#r|| -d s ,A] + /„ (M) + /j m) (M) , 



which are coupled by the constraint of inextensibility, 



i(d s rj 2 + 0(|d s rj 4 ). 



(4) 



For simplicity in most of this paper, we consider the 
Rouse model which assumes local friction. We also 
focus on dilute solutions. Long range hydrodynam- 
ics and non-zero concentration of filaments will modify 
some of the results (see below). For the Rouse model 
h{s - s') = S(s - s'VCi where Ci = 4^ = 2Q| - We 
model the active velocities by a Gaussian noise with mean 



<fj m) > = and (/. 



a drift, reflecting the po- 



lar nature of the filament and mean square fluctuations, 

Sf (m) = f (m) _ ( f (m)} giyen by 

(Sft l) (s,t)Sf { ; n \s',t')) = yffiQSaSis - s'Mt t>) 

(5) 

where = exp{— \t\/r} and {i,j} refer to _L, ||. The 
level of activity is controlled by the parameter O. Cor- 
relations of the active force decay over a time, r the 
typical activity period. The constants ctij measure the 
relative partitioning of the activity between transverse 
and longitudinal components and satisfy the relationship 
a 2 | + a\ = 1. In general, the force applied by the mo- 
tor protein on the filament will not be purely tangential 
depending for example on its trajectory of approach of 
the filament and the conformation of the neck and chain 
region [jjj. Averaging over orientations, only the longitu- 
dinal component will have a non-zero average, because of 



the filament polarity. In addition, the remainder of the 
energy of the actomyosin reaction that is not converted to 
directed work will be dissipated as heat contributing to 
the fluctuations on the same (reaction) time-scale. It is 
difficult apriori to estimate on^ and we choose somewhat 
arbitrarily, aj_ = aii = \j\[2. Between active events the 
'motors' diffuse freely in the solution so it is reasonable 
to assume that there is no spatial correlation between ac- 
tive sites. We emphasise that unlike the thermal noise, 
the active noise correlations do not satisfy the fluctuation 
dissipation theorem. 

The eqns.(||j3|^) are most easily studied by analysing 
the motion of bending modes of wave-vector q and fre- 
quency uj. Defining F(s, t) = f ^^F(q, u) exp(iut+qs) 
we only consider q such that -k / a q 3> tt/L On length- 
scales below L p , Ct(s) can be obtained from the trans- 
verse fluctuations. The transverse dynamics are approx- 
imately given by 



f ± ( g ,^+f| m) (g,^) 
iuj + aq 4 



(6) 



where a — k/Cl- Corrections due to the tension A(s,t) 
are higher order in the gradient expansion. The trans- 
verse fluctuations at time t are relaxed over a length- 
scale £j_(t) = (at) 1 / 4 . Because of inextensibility, eqn. 
(^) there is an induced time-independent longitudinal 
motion r^p (s) due to the averaged transverse motion, 

5 s rj 0) ( S ) = (i(9 s r ± ) 2 ) ± . 

The tangent correlation function is given by 



2 



Ct (s) = (t(s, t) ■ t(0, t)) ~ (l - 9 s rj 0) (a, i) - d s r j 0) (0, i) + S s r ± (s, i) • 5 s r ± (0, t)) 



(7) 



Using eqns. (f|,|]j7]), the tangent correlation function 
has the form 



C t (s) ~ exp 



s 

£11 



2Q T dq 1 — cos((7s) 
^"7 2^g 2 ( 9 4 + (ar)- 1 ) 



(8) 



The correlation function is plotted in Fig. J|. For length- 
scales less than £ c ~ = (kt/C^) 1 ^ the effective 
persistence length is approximately equal to the bare 
persistence length L p whilst for length-scales above £ c , 
the conformations can be well modelled by a worm- 
like chain with a lower renormalized persistence length 
given by L* = L p (1 + Qr^±/kBT)~ . Dimensional anal- 
ysis suggests a new active temperature scale given by 
kBT a = ®t(±. 

Dynamics The dynamics of inextensible filaments is 
anisotropic pTj| , |il^ , p3||T4| ] : longitudinal motion has differ- 
ent relaxation dynamics to the transverse. Viscous dis- 
sipation due to longitudinal motion must be taken into 
account, giving a time-scale over which 'tension' propa- 
gates along the filament 0|l]|J||r|. From eqn. (§), 
we can calculate the transverse dynamic fluctuations, 



R\(t) 



i2^3/4 C -3/4 r7/4j i<<r 

t 3/4£M,i/4 Ci -3/^ BT + q Ut)> t » r 



(9) 



There is a cross-over from ballistic motion at short times 
to sub-diffusive behaviour at long times. Note the ther- 
mal and active contributions to the sub-diffusive regime. 

We can solve for the longitudinal motion self- 
consistcntly as follows: we average over only the 
transverse fluctuations and use the inextensibility con- 
straint (l^Jl^l to give a relationship between the ten- 
sion and the longitudinal motion, iq(r\\ (q, u>) — (q)) — 

K(oj) A(g, to) defining a frequency dependent extensional 
compliance, K(lo) = K eq (uj) + K a (u) which has equilib- 
rium and active contributions. The equilibrium modulus 
is given by K eq (uj) ~ 2~^ 4 k B TK,- 5 / 4 (i(±ojy 3/4 is not 
new and has been obtained previously |ll|Jl^]. We find 
in addition a new active contribution given by 



K a (L0) 



49r r (ak i r + f(ak 4 r) - iujt) 

~ J k f{ak 4 r){2ak 4 - iu){f(ak 4 T) - icor) 

(10) 



where f(x) = x + 1. 



In the limit to — > 0, it is given by 

[^a)- 3/4 -i(a/2r)- 3 ^ 

O(uo). The modulus K(u>) is then substituted into the 
longitudinal dynamics. This self-consistent approach cor- 
responds to an infinite resummation of a set of diagrams 



of the perturbation expansion in A JL3J . We obtain the 
following equation for the longitudinal motion, 

- ( \ -Q) , x /||(g^) + /|{ m W) f . 

r\\{q,w -r '(q)~ , " (11) 

11 iuj + q z K 1 (w)/C|| 

As defined above, r^ (s) is the time-independent motion 
in the longitudinal direction due to the averaged trans- 
verse motion. Given this, we calculate the longitudinal 
dynamical fluctuations, -Rjj(i) = \(ri\(s,t) — rn (s, 0)) 2 )„ 



t 2 (yl + 6C|! /2 2- 3 / 4 K 3 / 8 ci 3/ V /8 ) , t 

f 2„,2 i ,7/8 8r(l/8) (2 K ) 5 /'(t s T+e( 1 r)^ 8 . 
b u m ' b 7tt ,7/8 l " 



< T 
T . 



(12) 



Because of the drift term (the directed motor/filament 
interaction), the longitudinal diffusion is ballistic for all 
time-scales. 

High frequency viscoelasticity The complex shear 
modulus of a solution of semiflexible filaments 
at high frequencies is dominated by the exten- 
sional compliance and G*(u>) — j^pK^ 1 ^) ~ 
0.133/V (K eq (uj) + K a (uj)) §,|n],|l2| which for a solution 
of passive filaments is given by G*(u>) oc ui 3 ^ 4 . For the 
active filaments, using eqn. (Q, we plot the absolute 
value of the complex high frequency modulus for vary- 
ing activity in Fig. |^. We see a crossover from at high 
frequencies a modulus corresponding to the bare tem- 
perature, T to at low frequencies a modulus equivalent 
to the higher renormalized active temperature T + T a , 
both with a scaling of uj 3 ^ 4 . In the cross-over regime, the 
modulus appears to have a stronger frequency depen- 
dence with a power law G*(lu) oc ui a where a > 3/4. It 
is interesting to note that the cross-over occurs over a 
very wide frequency range. 

Hydrodynamics Within the (screened) Oseen approx- 
imation the mobility tensor is given by |t]] , Hij [r] = 

%Jt,\t\ (ftj + ^ r i) for l r l > a - Tne hydrodynamic 
screening length £ — > oo for dilute solutions and is 
equal to the mesh-size £ = (p a a) -1 / 2 for semi-dilute 
solutions with actin concentration p a and filament di- 



ameter a.. Therefore, h(s — s' 
For dilute solutions (£ - 



— m eqns 



87r?7|s— s' 

oo), this gives loga- 
rithmic corrections to the Rouse model described above 
and implies a modification of the cross-over length to 
£ c ~ (kt/Q±_ log[(Kr/C_L) 1 ^ 4 /a]) 1 / 4 . For concentrated so- 
lutions (£ finite) , the friction coefficients cross-over length 
is 4~(WaiogK/a]) 1/4 - 
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Entangled solutions For entangled solutions, the 
filament can be modelled as confined in its tube. 
Given a mesh-size £, we can define a tube diame- 
ter D e ~ Lp(£/Lp) 6 / 5 and entanglement length L e ~ 
L p (£/L p ) 4 / 5 g|. The tube can be modelled as a con- 
fining potential for rj_ which we can model as Vmbe — 
^fc|rj_| 2 with k ~ chosen so as to give (\r±(L e ) — 

r^O)! 2 ) 1 / 2 = D e , the tube diameter. We obtain 
the tangent correlation function, Ct(s) = 1 — — 

2© f dq_ 1 — cos(qs) 

a 2 J 2tt q 2 {q 4 - + {aT)- 1 +k/n) ' 

Discussion Our results should be relevant for recent 
experiments on mixtures of F-actin and SI Myosin [Q. 
The molecular motor myosin interacts with F-actin in the 
presence of ATP undergoing a conformational change in 
the process. Hydrophobic interactions between tails of a 
common variant, Myosin II, lead to multi-headed clusters 
which can act as active cross-links [|]|| between two or 
more filaments. In contrast, SI myosin is a single-headed 
version of myosin that is without a tail. Therefore, in 
general SI myosin/ ATP interacts with single polar actin 
filaments applying biased non-equilibrium forces. The SI 
experiments found surprisingly different steady-state and 
dynamic behaviour of the filaments as compared to the 
pure F-actin system 

We now estimate the values for the parameters of our 
model corresponding to the SI experiment. For an actin 
monomer concentration p a , the fraction of bound myosin 
4> can be estimated using the equilibrium constant |Q, 
k eq — 500 nM for the passive reaction (Actin + Myosin 
Actomyosin) . The fraction of bound myosin is given by 
4> = Pal {k eq + p a ) leading to the typical separation of the 
active sites (motors) on the filament, i m ~ ((f'PmS, 2 )^ 1 = 
a(p a /4>Pm) where p m is the concentration of SI myosin. 
Then we may divide the filament into regions of size £ m 
in which the motors are expected to act independently. 
Therefore (f{s,t)f{s',t')) = £ m S(s - s')<f m (t)fm (*')>■ We 
assume that motor attachment on the filaments is a Pois- 
son process |]l7| , i.e. that motors arrive at random times 
t n with a constant rate A. The number of motors arriving 
in a period At ^> 1/A has a Poisson distribution. The 
forces applied by each motor are assumed to decay over 
a timescale r so that the force at time t is given by 

F m (t) = J2j Lf o9(t-t n ) (13) 

n— 1 

where g(t) = exp(— f/r), t > and g(t) = 0, t < and fo 
is the typical force applied by a motor to a filament of 
diameter a. After some standard manipulations JTfl] and 
averaging over orientations, the velocity correlations can 
be shown to be given by eqn. (|^) above with the activity 
parameter given by O ~ Arf" 1 (/oa) 2 and we set A = 1 /r. 
Similarly, the drift can be estimated as v m ~ (/o/C||)^m • 
The typical force, fo and active time-scale, r can be ob- 
tained from biochemical data. The viscosity of water is 
7] ~ 10~ 3 Pa s and F-actin has diameter a ~ 7nm and 



persistence length L p ~ 15/im. The stall force of Myosin 
is approximately 5 pN and the duration of the motor cycle 
is approximately r ~ 5 ms The motor step-size d m is ap- 
proximately d m ~ 10 ran. For ATP saturation and 5/iM 
SI and 14/xM F-actin concentrations, we find cf> w 0.9 and 
l m ss 2a. This leads to an estimate for 0£j_t ~ ksT at 
310 K giving L* ~ \L p and £ c ~ 0.1pm. 

In conclusion, we have studied a simple model of ac- 
tive filaments and obtained different static and dynami- 
cal properties as compared to passive semiflcxible poly- 
mers. These differences should be observable with video 
microscopy and in linear rheological experiments and be 
relevant for experiments on actin-myosin systems. 
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FIG. 1. (a) The fluctuating polar filaments of persistence 
length L p parametrised by R(s, t) decorated by active centres. 
Note that there is a (+) and a (— ) end for each filament, (b) 
The filament on length-scales below L p showing the transverse 
r±(s,t) and longitudinal r||(s,f) motion, (c) A schematic of 
the cycle of activity of a motor such as myosin with activity 
time t. 



FIG. 2. The tangent correlation function for the active 
filament of bare persistence length L p = 10, compared with 
those simple wormlike chains of persistence length of L v = 10 
and Lp — 5 using units such that ar = 1. Inset: Ct(s) near 
s = 0. 
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FIG. 3. The magnitude of the high frequency modulus for 
active filaments. We have chosen units of length and time such 
that ar = 1 and feT/d = 1. In these units, the bare persis- 
tence length of the filaments is L p = 1. (a) The cross-over to 
a renormalised persistence length Lp = 1/101. (b) Increasing 
activity corresponding to Q(±t = ksT a = 25, 50, 100 
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